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Abstract. We determine the limit distributions of sums of deterministic chaotic variables in 
unimodal maps assisted by a novel renormalization group (RG) framework associated to the 
operation of increment of summands and rescaling. In this framework the difference in control 
parameter from its value at the transition to chaos is the only relevant variable, the trivial fixed 
point is the Gaussian distribution and a nontrivial fixed point is a multifractal distribution with 
features similar to those of the Feigenbaum attractor. The crossover between the two fixed 
points is discussed and the flow toward the trivial fixed point is seen to consist of a sequence of 
chaotic band mergers. 



1. Introduction 

The hegemony of the Central Limit Theorem [TJ [2] for sums of deterministic variables 
generated by a number chaotic mappings have for some time been observed and also 
mathematically proved [3]. Since the mixing properties of chaotic trajectories yield variables 
indistinguishable to independent random variables, it is of interest to study nonmixing systems 
such as mappings at the transition from regular to chaotic behavior. Recent numerical 
explorations of time averages of iterates at the period-doubling transition to chaos [8] have been 
presented and interpreted as possible evidence for a novel type of stationary distribution. 

The dynamics toward and at the Feigenbaum attractor is now known in much detail [9j [10] , 
therefore, it appears feasible to analyze also the properties of sums of iterate positions for this 
classic nonlinear system with the same kind of analytic reasoning and numerical thoroughness. 
Here we present the results for sums of chronological positions of trajectories associated to 
quadratic unimodal maps. We consider the case of the sum of positions of trajectories inside 
the Feigenbaum attractor as well as those within the chaotic 2 n -band attractors obtained when 
the control parameter is shifted to values larger than that at the transition to chaos. From 
the information obtained we draw conclusions on the properties of the stationary distributions 
for these sums of variables. Our results, that reveal a multifractal stationary distribution 
that mirrors the features of the Feigenbaum attractor, can be easily extended to other critical 
attractor universality classes and other routes to chaos. 

The overall picture we obtain is effectively described within the framework of the 
renormalization group (RG) approach for systems with scale invariant states or attractors. 



Firstly, the RG transformation for the distribution of a sum of variables is naturally given 
by the change due to the increment of summands followed by a suitable restoring operation. 
Second, the limit distributions can be identified as fixed points reached according to whether 
the acting relevant variables are set to zero or not. Lastly, the universality class of the non-trivial 
fixed-point distribution can be assessed in terms of the existing set of irrelevant variables. 

As it is well known [8] a few decades ago the RG approach was successfully applied to the 
period-doubling route to chaos displayed by unimodal maps. In that case the RG transformation 
is functional composition of the mapping and its effect re-enacts the growth of the period 
doubling cascade. In our case the RG transformation is the growth and adjustment of the 
sum of positions and its effect is instead to go over again the merging of bands in the chaotic 
region. 

Specifically, we consider the Feigenbaum map g(x), obtained from the fixed point equation 
g{x) = ag{g{x/a)) with g(0) = 1 and g'(0) = 0, and where a = —2.50290... is one of 
Feigenbaum's universal constants [8j. For expediency we shall from now on denote the absolute 
value \a\ by a. Numerically, the properties of g{x) can be conveniently obtained from the 
logistic map /^{x) = 1 — fix 2 , — 1 < x < 1, with fi = = 1.401155189092... The dynamics 
associated to the Feigenbaum map is determined by its multifractal attractor. For a recent 
detailed description of these properties see [9[ HO] . For values of fi > [i^ we employ a well- 
known scaling relation supported by numerical results. 

Initially we present properties of the sum of the absolute values \xt\ of positions xt = 
ffj,oa,2( x t-i)i t = 1)2,3,..., as a function of total time N visited by the trajectory with initial 
position xq = 0, and obtain a patterned linear growth with N. We analyze this intricate 
fluctuating pattern, confined within a band of finite width, by eliminating the overall linear 
increment and find that the resulting stationary arrangement exhibits features inherited from 
the multifractal structure of the attractor. We derive an analytical expression for the sum that 
corroborates the numerical results and provide an understanding of its properties. Next, we 
consider the straight sum of xt, where the signs taken by positions lessen the growth of its 
value as N increases and the results are consistently similar to those for the sum of \xt\, i.e. 
linear growth of a fixed- width band within which the sum displays a fluctuating arrangement. 
Further details for the sum of xt are not included because of repetitiveness. Then, we show 
numerical results for the sum of iterated positions obtained when the control parameter is 
shifted into the region of chaotic bands. In all of these cases the distributions evolve after a 
characteristic crossover towards a Gaussian form. Finally, we rationalize our results in terms of 
an RG framework in which the action of the Central Limit Theorem plays a fundamental role. 

The starting point of our study is evaluation of 



with [i = fioo and with xo = 0. Fig. 1A shows the result, where it can be observed that the values 
recorded, besides a repeating fluctuating pattern within a narrow band, increase linearly on the 
whole. The measured slope of the linear growth is c = 0.56245... Fig. IB shows an enlargement 
of the band, where some detail of the complex pattern of values of y/M 00 (N) is observed. A 
stationary view of the mentioned pattern is shown in Fig. 1C, where we plot 
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in logarithmic scales. There, we observe that the values of y| ioo (A r ) fall within horizontal bands 
interspersed by gaps, revealing a fractal or multifractal set layout. The top (zeroth) band 




Figure 1. A) Sum of absolute values of visited points xt, t = 0,...,N, of the Feigenbaum's 
attractor with initial condition xq = 0. B) A closer look of the path of the sum (see dotted circle 
in A), for values of N around 67410. C) Centered sum y'(N) in logarithmic scales. See text. 



contains y'„ for all the odd values of N, the 1st band next to the top band contains y' for 
the even values of N of the form N = 2 + 4m, m = 0, 1, 2, ... The 2nd band next to the top band 
contains y| ioo (A r ) for iV = 2 2 + 2 3 m, m = 0, 1,2, and so on. In general, the fc-th band next to 
the top band contains y' flaa (2 k + 2 fc+1 m), m = 0, 1, 2, ... Another important feature in this figure 
is that the y' floo (N) for subsequences of N each of the form N = (21 + l)2 k , k = 0, 1, 2, with 
I fixed at a given value of I = 0, 1, 2, appear aligned with a uniform slope s = —1.323... The 
parallel lines formed by these subsequences imply the power law y'^ (JV) ~ N s for N belonging 
to such a subsequence. 

It is known j9j [TT] that these two characteristics of y' (N) are also present in the layout 
of the absolute value of the individual positions \xt\, t = 1,2, 3, ... of the trajectory initiated at 
xo = 0; and this layout corresponds to the multifractal geometric configuration of the points of 
the Feigenbaum's attractor, see Fig. 1 in [TT]. In this case, the horizontal bands of positions 
separated by equally-sized gaps are related to the period-doubling 'diameters' |S] set construction 
of the multifractal [10]. The identical slope shown in the logarithmic scales by all the position 
subsequences \xt\, t = (21 + l)2 fc , k = 0, 1,2, each formed by a fixed value of I = 0, 1,2, 
implies the power law \xt\ ~ t s , s = — In a/ In 2 = —1.3236..., as the \xt\ can be expressed as 
\ x t\ — \ x 2i+i\ a~ k ,t = (21 + l)2 fc , k = 0, 1, 2, or, equivalently, \xt\ ~ t s . Notice that the index 
k also labels the order of the bands from top to bottom. The power law behavior involving the 
universal constant a of the subsequence positions reflect the approach of points in the attractor 
toward its most sparse region at x = from its most compact region, as the positions at odd 
times |x2z+i| = x 2i+i, those in the top band, correspond to the densest region of the set. 

Having uncovered the through manifestation of the multifractal structure of the attractor into 
the sum y' (N) we proceed to derive this property and corroborate the numerical evidence. 
Consider Eq. @ with N = 2 k , k = 0, 1,2, the special case I = in the discussion above. 
Then the numbers of terms \xt\ per band in y^ ao (2 k ) are: 2 k ~ 1 in the top band (j = 0), 2 k ~ 2 
in the next band (j = 1),..., 2° in the (k — l)-th band, plus an additional position in the A;-th 



band. If we introduce the average of the positions on the top band 



(a) = 2-(*-D J2 xjy+i, (3) 

j=0 

the sum y l _ loo (2 k ) can be written as 

y, x (2 k ) = (a) 2*" 1 J2(2a)^ + a^ 1 ) + a k . (4) 

3=0 

Doing the geometric sum above and expressing the result as yu, ao (2 k ) = c2 k + da~ k , we have 

(a)a <a)2a^ + i _ 



2a -1 V 2a -1 

Evaluation of Eq. ^ yields to (a) = 0.8999..., and from this we obtain c = 0.56227... 
and d = 0.68826... We therefore find that the value of the slope c in Fig. 1A is properly 
reproduced by our calculation. Also, since In \y^ 00 (2 k ) -c2 k \ = lnd-kln a, or, equivalently, 

Iny'^N) = hid- JVlna/ln2, N = 2 k , k = 0,1,2,..., we corroborate that the value of the 
slope s in inset of Fig. 1C is indeed given by s = — In a/ In 2 = 1.3236... (We have made use of 
the identity a~ k = N~ lna / ln2 , N = 2 k , k = 0,1,2,...). 

We note that the sum of xt from t = to N = 2 k , i.e. considering the signs taken by positions, 
can be immediately obtained from the above by replacing a' 3 by (— lpa - - 7 as the xt of different 
signs of the trajectory starting at xq = fall into separate alternating bands (described above 
and shown in Fig. 1 of [II]). In short, x t ~ (-1)^+1 a~ j , t = (21 + l)2 k , k = 0, 1,2, ... As 
stated, our numerical and analytical results are in agreement also in this case. 

We turn now to study the sum of positions of trajectories when = /i— fj,^ > 0. We 
recall that in this case the attractors are made up of 2 K , K = 1,2,3, bands and that their 
trajectories consist of an interband periodic motion of period 2 K and an intraband chaotic 
motion. We evaluated numerically the sums y^N) for an ensemble of initial conditions xq 
uniformly distributed only within the chaotic bands, for different values of A//; y'^(N) was then 
obtained similarly to Eq. Q by substracting the average (yp(N)) X0 anci rescaling with a factor 
jy-i/2_ fp^g p ane i s i n Fig. 2 show the evolution of the distributions for increasing number 
of summands N for a value of A// (chosen for visual clarity) when the attractor consists of 
2 3 chaotic bands. Initially the distributions are multimodal with disconnected domains, but 
as iV increases we observe merging of bands and development of a single-domain bell-shaped 
distribution that as N — > oo converges in all cases to a Gaussian distribution. 

These numerical results can be understood as follows. We recollect [8] that the relationship 
between the number 2 K , K 3> 1, of bands of a chaotic attractor and the control parameter 
distance Afj, at which it is located is given by 2 K ~ A/x" K , K = ln2/ln<5 F , where S F = 0.46692... 
is the universal constant that measures both the rate of convergence of the values of \i at period 
doublings or at band splittings to [i^. For A^ small and fixed, the sum of sequential positions 
of the trajectory initiated at xq = 0, Eq. ([T]), exhibits two different growth regimes as the total 
time iV increases. In the first one, when N <C 2 K , the difference in value Sxt = xt(n) — xt(Hoo) 
between the positions at time t for [i and /Xqo do not affect qualitatively the multifractal structure 
of the sum at //oo, nor its associated distribution. This is because the fine structure of the 
Feigenbaum attractor is not suppressed by the fluctuations Sxt, as these contribute to the sum 
individually during the first cycle of the interband periodic motion. The discrete multi-scale 
nature of the distribution for fj,^ is preserved when the interband motion governs the sum 
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Figure 2. Distributions for the sums of \xt\, t = 0, ...,N, of an ensemble of trajectories with 
initial conditions within the 2 3 -band attractor at A/i = 0.0028448109. The number of summands 
N are indicated in each panel. See text. 

y^{N). In the second regime, when N 3> 2 K , the situation is opposite, after many interband 
cycles the fluctuations 5x t add up in the sum and progressively wipe up the fine structure of 
the Feigenbaum attractor, leading to merging of bands and to the dominance of the fluctuating 
intraband motion. Ultimately, as N — ► oo the evolution of the distribution is similar to the 
action of the Central Limit Theorem and leads to a Gaussian stationary end result. It is also 
evident that as A/j increases the first regime is shortened at the expense of the second, whereas 
when A/i, — > the converse is the case. Therefore there exists an unambiguous A/i-dependent 
crossover behavior between the two radically different types of stationary distributions. This 
crossover is set out when the 5xt fluctuations begin removing the band structure in y'^{N) when 
Afj, is small and ends when these fluctuations have broadened and merged all the chaotic bands 
and y'AN) forms a single continuous interval. When fx = fj,^ this process never takes place. 

We are in a position now to put together the numerical and analytical information presented 
above into the general framework of the RG approach. As known, this method was designed to 
characterize families of systems containing amongst their many individual states (or in this case 
attractors) a few exceptional ones with scale invariant properties and common to all systems 
in the family. We recall [12] that in the language of a minimal RG scheme there are two fixed 
points, each of which can be reached by the repeated application of a suitable transformation 
of the system's main defining property. One of the fixed points, is termed trivial and is reached 
via the RG transformation for almost all initial settings, i.e. for all systems in the family when 
at least one of a small set of variables, named relevant variables, is nonzero. To reach the other 
fixed point, termed nontrivial, it is necessary that the relevant variables are all set to zero, and 
this implies a severely restricted set of initial settings that ensure such critical RG paths. The 
nontrivial fixed point embodies the scale invariant properties of the exceptional state that occurs 
in each system in the family and defines a universality class, while the differences amongst the 



individual systems are distinguished through a large set of so-called irrelevant variables. The 
variables in the latter set gradually vanish as the RG transformation is applied to a system 
that evolves toward the nontrivial fixed point. Further, when any system in the family is 
given a nonzero but sufficiently small value to (one or more of) the relevant variables, the RG 
transformation converts behavior similar to that of the nontrivial fixed point into that resembling 
the trivial fixed point through a well-defined crossover phenomenon. The recognition of the RG 
framework in the properties of the sums of positions of trajectories in unimodal maps and their 
associated distributions is straightforward. It can be concluded right away that in this problem 
(as defined here) there is only one relevant variable, the control parameter difference A/x. There is 
an infinite number of irrelevant variables, those that specify the differences between all possible 
unimodal maps (with quadratic maximum) and the Feigenbaum map g{x). There are two 
fixed-point distributions, the trivial continuum-space Gaussian distribution and the nontrivial 
discrete-space multifractal distribution (as observed in Fig. 1C). As explained above, there is a 
distinct crossover link between the two fixed-point distributions. The RG transformation consists 
of the increment of one or more summands in the sum ([!]) followed by centering like in Eq. ^ . 
The effect of the transformation in the distribution of the sum is then recorded. For sums of 
independent variables the transformation is equivalent to the convolution of distributions. Our 
results correspond to the dynamics inside the atractors, however, if the interest lies in considering 
only the stationary distribution of sums that do not contain the transient behavior of trajectories 
in their way to the attractor |J] our results are expected to give the correct answers for this case. 

In summary, we have found that the stationary distribution of the sum of iterate positions 
within the Feigenbaum attractor has a multifractal structure stamped by that of the initial 
multifractal set, while that involving sums of positions within the attractors composed of 2 K 
chaotic bands is the Gaussian distribution. We have also shown that the entire problem can 
be couched in the language of the RG formalism in a way that makes clear the identification 
of the existing stationary distributions and the manner in which they are reached. These basic 
features suggest a degree of universality, limited to the critical attractor under consideration, in 
the properties of sums of deterministic variables at the transitions to chaos. Namely, the sums of 
positions of memory-retaining trajectories evolving under a vanishing Lyapunov exponent appear 
to preserve the particular features of the multifractal critical attractor under examination. Thus 
we expect that varying the degree of nonlinearity of a unimodal map would affect the scaling 
properties of time averages of trajectory positions at the period doubling transition to chaos, 
or alternatively, that the consideration of a different route to chaos, such as the quasiperiodic 
route, would lead to different scaling properties of comparable time averages. For instance, the 
known dependence of the universal constant a on the degree of nonlinearity z of a unimodal map 
would show as a z-dependent exponent s = In a/ In 2 that controls the scale invariant property 
of the sum of trajectory positions with xq = (shown in Fig. 1C ). 

We analyzed the nature and the conditions under which a stationary distribution with 
universal properties (in the Renormalization Group sense) occurs for sums of deterministic 
variables at the transition between regular and chaotic behavior, such as those studied here for 
dynamics at zero Lyapunov exponent. The nonexistence of fluctuations in such critical attractors 
implies a distribution of the sums of these variables strictly defined on a discrete multifractal set 
and therefore different from known (Gaussian or otherwise) continuum-space limit distributions 
for real number random variables. 
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